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ABSTRACT 

We  consider  the  spatial  Fourier  transform  for  wave  vector 

k  of  the  charge  distribution  of  the  electrons  in  a  plasma  with  particle 
density  *1 1  electron  and  ion  temperatures  T  and  T^>and  Debye  length 
D.  We  assume  the  absence  of  a  magnetic  field,  neglect  collisions,  and 
assume  nD^  »  1  .  The  statistical  average  of  calculated  as 

a  function  of  a  =  l/kD  assuming  complete  thermodynamic  equilibrium; 
that  component  of  |pi^gP»  which  keeps  in  phase  with  the  ion  charge  density 
fluctuations  is  also  calculated. 

The  frequency  spectrum  of  the  time -varying  function  is  ob¬ 

tained  at  thermal  equilibrium  and  simplified,  assuming  the  ion  mass  to 
be  much  larger  than  the  electron  mass,  for  general  values  of  a  and 
T/Tp  For  small  a  the  main  component  of  the  spectrum  has  the  char¬ 
acteristic  Doppler  broadening  shape  corresponding  to  the  electron's  thermal 
velocity.  For  large  a  we  have  a  component  with  narrow  width  corres¬ 
ponding  roughly  to  the  ion-velocity  Doppler  spread  and  very  narrow  side 
bands  at  plus  and  minus  the  frequency  of  electrostatic  plasma  oscillations.* 

I.  INTRODUCTION 

In  the  last  decade  or  two  many  calculations  have  been  carried  out 
on  the  time  development  of  fluctuations  of  charge  density  in  an  ionized  gas 
under  a  variety  of  conditions.  For  a  given  volume  V  containing  N  electrons 

and  pj^.  have 

been  introduced  (mainly  for  mathematical  convenience),  which  are  essentially 

*This  paper  is  to  be  published  in  the  Physical  Review.  Derivations  in  Appendix 
I  and  II  were  prepared  by  P.  Goldreich. 


and  N/Z  positive  ions  of  atomic  charge  Z,  quantities  pj^^ 
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the  spatial  Fourier  transforms  for  wave  vector  k  of  the  electron  and  ion 
charge  densities.  More  specifically 


where  r>,(t)  and  ^.(t)  are  the  positions,  as  a  function  of  time  t,  of  the 
“J  J 

th 

j  electron  and  ion,  respectively,  and  k  is  a  constant  wave  vector. 

The  use  of  backs cattering  of  a  radar  beam  from  the  ionosphere  at 
great  heights  or  from  the  exosphere  has  been  proposed  recently  by  Gordon^ 
for  measuring  electron  density  and  temperature  at  various  heights.  Radar 
frequencies  of  50  to  1000  Mc/s  are  used  which  are  very  large  compared 
with  the  electron  plasma  frequency  (of  the  order  of  0.  1  to  10  Mc/s), 

Wp  =  (4irne^/m)^  ,  (2) 

where  n  =  N/V  is  the  particle  density  of  the  electrons  alone,  e  is  the 
charge  (in  C.  G.  S.  units)  and  m  the  mass  of  the  electron.  Standard  magne- 
toionic  theory  replaces  the  electrons  by  a  continuous  medium,  whose  re¬ 
fractive  index  is  close  to  unity  at  these  high  frequencies  and  would  not 
lead  to  reflection  or  attenuation,  if  the  density  of  the  medium  is  assur  led 
to  be  smooth  and  varying  slowly.  The  actual  amplitude  of  radiation  scattered 
for  a  scattering  angle  6  from  a  volume  V  of  ionized  gas  (dimensions  of  the 
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order  of  1  km)  is  then  simply  the  sum  of  the  Thomson  scattering  amplitudes  • 
from  each  of  the  N  electrons  in  the  volume  (Thomson  scattering  from  the 
positive  ions  is  negligible  because  of  their  large  mass).  Since  the  electrons 
are  highly  nonrelativistic,we  can  neglect  retardation  effects  and  the  scattering 
amplitude  from  each  electron  contains  a  phase  factor  like  those  appearing  in 
Equation  (1)  with  the  wave  number  k  given  by  k  -  4TTsin  j0/\  where  X. 
is  the  wavelength  of  the  electromagnetic  radiation.  The  total  amplitude  of 
backs catter  as  a  function  of  time  is  thus  proportional  to  the  quantity  g(t)  , 
defined  in  Equatioi^' (1),  for  a  fixed  value  of  k.  For  low  enough  intensity 
of  the  radar  beam  and  for  radar  frequency  large  compared  with  Wp  we  can 
neglect  altogether  the  effect  of  the  electromagnetic  radiation  on  the  quantity 

‘^k.  e 

Let  D  be  the  Debye  length,  defined  for  the  electrons  alone;  then 

2  - 

D  =  (,<T/4une  )*  ,  (3) 


where  n  is  the  electron  particle  density  and  T  the  electron  temperature, 
and  let  A  be  the  dimensionless  ratio 

A=nD^  oc  (e^n^/^/ZCT)"^^^  oc  D/CX/e^  .  (4) 

We  shall  only  consider  cases  throughout  this  paper  where  A  »  1  ;  i.  e.  , 

where  a  sphere  of  radius  equal  to  the  Debye  length  contains  very  many 

electrons  and  where  the  Coulomb  interaction  energy  between  "nearby" 

- 1  /3 

electrons  (separations  ~  n  '  )  is  small  compared  with  the  thermal  energy 
KT  .  This  inequality  certainly  holds  for  densities  and  temperatures  encountered 
in  the  ionosphere  and  exosphere,  as  well  as  for  many  laboratory  experiments. 
We  further  assume  throughout  that  the  gas  is  ionized  enough  so  that  collisions 
with  neutral  gas  atoms  or  molecules  can  be  neglected.  This  assumption  is 
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not  valid  for  the  lower  ionosphere  but  holds  for  heights  of  about  300  km  and 
higher.  For  a  highly  ionized  gas  with  A»  1  the  main  collision  process 
for  electrons  and  ions  is  multiple  Coulomb  scattering  through  small  angles, 
and  the  effective  mean  free  path  I  for  appreciablfs;  deflections^  is  of  the 
order  of  DA/logA»  D.  We  shall  assume  that  A  is  sufficiently  large  for 
the  mean  free  path  I  to  be  large  compared  also  with  our  effective  scale 
length  k“^  and  we  shall  neglect  collisions  altogether.  We  shall  calculate 
in  this  paper  only  in  the  absence  of  any  magnetic  field.  For  the 
ionospheric  applications  the  neglect  of  the  earth's  magnetic  field  is  not 
justified  at  the  lower  frequencies  of  about  50  Mc/s  and  is  expected  to  give 
a  moderately  good  approximation  at  radar  frequencies  of  about  400  Mc/s  or 
higher.  In  many  calculations  in  this  paper  we  assume  complete  thermo¬ 
dynamic  equilibrium,  but  we  shall  also  discuss  some  limited  deviatidna.^ 
from  such  equilibrium. 

In  all  our  calculations  the  following  dimensionless  parameter  a 
will  be  of  importance 


==  ^  ==  (4irne^/k^^T)^  =  X.  • 


kD 


W 


(5) 


In  the  limit  of  o  -►  0  the  collective  effects  of  the  Coulomb  interactions, 

which  become  important  only  over  distances  as  large  as  the  Debye  length 

D  or  larger,  are  negligible  over  distances  as  small  as  the  scale  length  k~^  , 

and  the  electrons  are  randomly  distributed  in  space.  In  this  case  we  have 

completely  incoherent  scattering  from  each  of  the  electrons,  at  least  at 

thermodynamic  equilibrium,  i.  e.  ,  in  evaluating  >  all  the  cross- 

2 

terms  average  to  zero;  and  we  simply  obtain  Ne  for  this  quantity.  In 
the  ionospheric  applications  for  a  radar  frequency  of  450  Mc/s,  for  in-  . 
stance,  a  <  1  at  heights  of  1000  or  2000  km  and  higher,  and  this  case  of 
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a  0  has  been  treated  in  detail  by  Gordon.  ^  We  shall  carry  out  calculations 
in  this  paper  for  arbitrary  values  of  a  for  which  the  spatial  correlations 
between  all  the  electrons  and  ions  have  to  be  taken  into  account. 

In  Section  II  we  evaluate  the  time  average  of  the  intensity  of  the 


electron  density  fluctuation,  jPi^e  |  »  arbitrary  values  of  the  para¬ 

meter  a  and  the  atomic  charge  Z  of  the  positive  ions  at  complete  thermo¬ 
dynamic  equilibrium.  We  shall  calculate  these  averages  from  first  principles, 
although  they  could  be  obtained  more  easily  using  results  from  the  Debye - 
Huckel  theory.  Such  a  calculation  does  not  give  the  time  development  of 
Pke(t)  or  its  frequency  Fourier  transform,  but  the  intensity  can  be  divided 
into  two  parts  with  different  characteristic  frequency  spreads,  if  the  ion 
mass  M  is 'very  large  compared  with  the  electron  mass  m  .  The  ions  move 
very  slowly  compared  with  the  electrons,  and  if  we  consider  the  ions  fixed, 
we  can  evaluate  the  average  correlation  of  the  electron  density  distribution 
with  that  of  the  ions.  Such  a  calculation  will  give  that  part  of  that  vpries 
very  slowly  with  time  (characteristic  of  ion  thermal  velocities).  The  re¬ 
maining  part  of  pj^g  varies  rapidly  with  time  (characteristic  of  electron 
thermal  velocities). 

In  Section  III  we  derive  formulas  for  ^|Qj^^(w)p^  where  Qj^g(w)  is 
the  frequency  Fourier  transform  of  limit  of  7  -*■  0  .  This 

quantity  is  relevant  if  a  frequency  spectrum  is  observed  over  a  long  but 
finite  time  period.  The  calculations  are  carried  out  for  complete  thermo¬ 
dynamic  equilibrium  except  that  the  electron  and  ion  temperatures  T  and 
Tj^  need  not  be  equal.  In  Section  IV  the  general  results  are  simplified  and 
approximations  evaluated,  using  the  fact  that  the  tion  mass  M  is  much 
larger  than  the  electron  mass  m  .  In  Section  V  a  special  kind  of  deviation 
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from  equilibrium  is  discussed,  where  we  assume  that  an  external  agent 
suddenly  alters  the  degree  of  ionization  in  a  nonuniform  manner  at  some 
time,  but  the  medium  is  allowed  to  relax  to  equilibrium  after  this  time. 

II.  SOME  TIME -AVERAGED  INTENSITIES 


We  consider  a  volume  V  containing  N  electrons  and  N/Z  positive 
ions  of  atomic  charge  Z  with  N  and  V  extremely  large  but  with  electron 
particle  density  n  =  N/V  fixed  and  finite.  In  this  section  we  assume  comr 
plete  thermodynamic  equilibrium  at  temperature  T  and  evaluate  statistical 
averages  of  various  quantities.  Here  we  need  make  no  assumptions  about 
the  collision  mean  free  path  being  large  as  we  do  in  the  remaining  sections. 
We  do  have  to  assume,  however,  that  the  dimensionless  parameter  A  defined 
in  Equation  (4)  is  large  compared  with  unity;  in  an  expansion  in  inverse 
powers  of  A  we  shall  calculate  explicitly  only  the  leading  term  and  give  only 
qualitative  estimates  of  higher  order  corrections. 

Let  k  be  any  wave  vector  which  satisfies  periodic  boundary  con¬ 
ditions  for  the  volume  V  and,  for  any  given  spatial  distribution  of  all  the 
electrons  and  ions,  define  complex  quantities  ,  pj^-  and  pj^^  according 

to  Equation  (1).  Except  for  the  relation  p_j^  =  .  different  values  of  k 

represent  independent  modes.  We  shall  write 


-i6 


^ke  ^'’ke® 


ke 


Pki  =  Pki® 


ki 


-i£ 


Pkt  =  Pke  +  Pki  =  Pkt« 


kt 


where  each  p,  and  6  is  real  and  positive.  We  shall  need  the  electrostatic 

potential  <t>(r)  and  electric  field  E(r)  arising  from  the  given  distribution  of 

electrons  and  ions.  After  carrying  out  a  Fourier  transformation  of  the 

Coulomb  potential  -e/|'r  -  r.  resulting  from  the  electron  (and  of  its 

3 
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gradient)  and  summing  over  all  electrons  and  ions,  we  find 


COS  (k»r  -  6  )  , 

k 

E(£)  =^k  (Bir/Vk^)  jij^^sin(k*£  -  6^^)  ,  (7) 

k 


where  the  k  summation  is  carried  only  over  half  of  all  the  possible  k 
vectors  (those  with  positive  z  component,  say). 

We  now  consider  N-1  of  the  electrons  and  all  the  ions  as  fixed, 
introduce  an  electron  and  ask  for  the  statistical  ensemble  average  of 

cos(k»r  -  6)  over  all  positions  r  of  this  extra  electron,  k  and  6  being 
fixed.  This  average. is 


f  3  e4.(r)/AC  T 

J  d'^r  cosj^  r^  6)e  _ 


<cos(k.r  -  6)>  =  —  —'e^JrjyirT' 
'  d  re 


I 

^  ■ 


(8) 


where  K  is  Boltzmann's  constant.  As  V  -►  oo  we  could  write,  no  matter 
what  the  value  of  A  , 
e<t>(r)//<  T 


=UD  +  (Sire/Vk^/C  T)^i^^cos(k.r  -  . 


(9) 


and  rewrite  this  infinite  product  as  an  infinite  series  of  terms  with  successive 
positive  powers  of  l/ KT  .  However,  this  seties  will  converge  rapidly  if, 
and  only  if,  je(f)(£)//C  T  |«  1.  If  A  »  1  the  Coulomb  interaction  between 
"neighboring"  particles  is  weak  compared  with  /^T  and  this  inequality  is 
satisfied  for  all  values  of  r  except  those  very  close  to  one  of  the  fixed 
charges.  In  the  integral  in  Equation  (8)  distances  away  from  a  fixed 
pharge  that  are  small  compared  with  both  the  scale  length  k“^  and  the 


-7- 


Debye  length  D  are  unimportant.  At  a  distance  k"^  ,  the  potential  energy 
of  an  electron  is  of  order  e  k  and  the  use  of  an  expansion  in  powers  of 
1//CT  will  give  rapid  convergence  as  long  as  (Aa)"^  ~  e^k//CT  «  1  ,  as 
well  as  A  »  1  .  We  shall  assume  that  this  inequality  also  holds  (it  breaks 
down  only  for  a  «  A  1  ,  and  for  such  very  small  values  of  a  the 
Coulomb  correlations  are  negligible  and  the  problem  trivial  in  any  case). 
Using  the  expansion  of  Equation  (9)  in  the  integrands  of  Equation  (8),  we 
keep  only  the  zero  order  term  (unity)  in  the  denominator.  In  the  numer¬ 
ator  the  zero  order  term  gives  no  contribution,  we  keep  only  the  terms  of 
first  order  in  l/K.'T  ,  and  since  the  cosine  terms  for  different  values  of 
k  are  orthogonal,  we  obtain 

<cos(krr  -  6)>  =  (4ire/Vk^/<T)ji  cos(6  -  6,  J  ,  (10) 

—  —  kt 

where  the  and  refer  to  N-1  electrons  and  N/Z  ions.  A 

similar  calculation  for  <cos  (k*R  -  6)  >  ,  where  R  is  the  position  of  an 
additional  positive  ion,  simply  gives  -Z  times  the  expression  in  Equation 


The  use  of  Equation  (1)  gives  a  double  sum  over  indices  j,  I  for 
2 

a  quantity  like  for  all  N  electrons.  In  this  double  sum 

we  separate  out  the  terms  with  j  =  Z  for  which  the  phase  factors  cancel. 

When  a  statistical  ensemble  average  is  taken,  the  various  terms  in  the 

remaining  double  summations  give  identical  results,  and  replacing  N(N-l) 
2 

by  N  and  dropping  the  subscript  k,  we  obtain 


<p.^>  =  Ne^l^l  +  N<cos[jc*(rj  -  £2)^]^^  » 

<M^>  =  Ne^l  Z  +  N<cos[;;k.(Rj  -  R^)l>|  . 

<jx^>  r  <ji2>  +  <^2>  _  2N^e^  <cosk.(r.  -  R^)>  , 

J 
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where  and  R  j  are  the  positions  of  the  electron  and  ion  respectively. 

The  cosine  expectation  values  are  the  Fourier  transforms  of  the 
two -particle  correlation  functions,  which  would  vanish  for  randomly  dis¬ 
tributed  electrons  and  ions,  and  we  evaluate  them  as  follows.  We  keep 
for  all  particles  except  electrons  j  and  I  fixed,  keep  r  .  fixed  at 
first,  and  average  over  r  ^  .  This  average  can  be  obtained  in  analogy  with 
the  derivation  of  Equation  (10)  but  in  Equation  (9)  we  have  to  add  the  term 

Ql  -  (8ire^/Vk^/CT)  cos  k*(£  •  -  £/)I]  • 

which  represents  the  Coulomb  interactions  between  the  two  electrons  singled 
out.  We  have 


<cos  k*(£j  -  £^)>^  =  (4ire/Vk  /C  T)  Q-e  +  cos  (k*£j  -  6^)3]  . 

J 

with  r  .  fixed.  We  average  next  over  r  .  ,  still  keeping  p  fixed,  use 
“  J  “  J  t 

Equation  (10),  and  finally  average  over  the  remaining  particles  as  well. 
Using  the  same  procedure  on  the  other  cosine  terms  we  find 

-2  -1 

<cos  k*(r.  -  r,)>=  Z  <cos  k»(R  .  -  R  )>  =  -Z  <cos  k*  (r  .  -  R  ,  )> 

-  -J  -4  -  J  -  '-J  —I' 

=  (aVN)Ca^<|ij>^  (Ne^)-‘  -  ij  .  (12) 

2 

where  the  diniensionless  parameter  a  is  defined  in  Equation  (5); 
refers  to  N-2  electrons,  but  N-2  may  be  replaced  by  N  ,  since  we 
expect  no  cancellation  of  large  terms. 

2 

Substituting  Equation  (12)  into  Equation  (11)  to  eliminate 

2  2 

and  <|i.  >  yields  an  explicit  expression  for  <}ji.  >  ,  and  substituting  this 

I  »  , 

expression  back  into  Equation  (12)  gives  explicit  expressions  for  the 
cosine  expectation  values.  Using  Equation  (11)  again,  we  finally  obtain 
the  desired  expressions; 


-9- 


(13) 


<|ig>  =  Ne^  (1  +  Za^)D  +  (Z  +  l)a^T^ 

<|i?>  =  ZNe^(l  +  a^)Dl  +  (Z  +  l)a^T^ 

<\i^>  =  (Z  +  l)Ne^D  +  (Z  +  l)a^T^ 

N<cosk- (r  j  -  r  ^)>  =  -a^Cl  +  (Z  •»•  l)a^T^  •  (14) 

Equation  (14)  merely  represents  a  rederivation  of  the  Fourier  transform 
of  the  well-known  Debye -Hiickel  two -particle  correlation  function  for  two 
electrons.  If  we  had  taken  over  this  expression  (and  similar  ones)  from 
the  Debye -Hiickel  theory,  substitution  into  Equation  (11)  would  have  given 
the  desired  results  without  requiring  any  other  formulas  of  the  present 
section. 

For  a  «  1  our  results  in  Equation  (13)  reduce  to  those  for  ran¬ 
domly  distributed  particles,  ad  they  should.  For  a  «  1  the  total  charge 

2 

density  fluctuations  <\i^  >  are  smaller  than  those  for  random  distribution 
by  a  factor  a  cck  .  The  electron  charge  density  fluctuations  » 

however,  are  reduced  only  by  a  factor  Z/(Z  +  1)  even  in  the  limit  of 

•  3  2 

a  00  .  Pines  and  Bohm  have  carried  out  calculations  for  <^g>  for  a 

model  in  which  the  positive  charges  are  uniformly  and  continuously  dis¬ 
tributed.  Their  results  can  be  obtained  from  our  more  general  ones  in 

Equation  (13)  by  making  the  formal  substitution  Z  =  0  ,  in  which  case 

2  2  2  2 
<|Ag>  =  <|i^  >  and  <|i.  >  =  0  .  The  expression  for  ^  Z  =  1  , 

the  case  of  greatest  interest  for  the  ionospheric  application,  has  also  been 

derived  by  different  methods  by  Fejer^  and  by  Renau.  ^  Kahn^  has  derived 

2 

an  expression  for  general  values  of  Z  and  a,  which  agrees  with 

ours  in  the  two  limiting  cases  a  «  1  and  a  »  l.but  appears  to  be  in- 
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correct*  for  general  values  of  a  . 

We  finally  evaluate  in  a  similar  manner  another  statistical  average 
which  forms  only  part  of  •  We  first  consider  the  positions  of  the  ions, 

and  hence  p.  =  |i.e  ,  as  fixed  and  ask  for  the  average  of  the  component 
of  the  electron  quantity  p  which  is  in  phase  with  the  constant  p.  .  We 
first  rewrite  the  expression  in  Equation  (10)  for  the  average  over  the 
electron  (with  N-1  electrons,  as  well  as  the  ions,  fixed),  using  the 
definitions  in  Equation  (6): 

<cos  (k^r  -  6)  >  =  (4Tre/Vk^/<[T)Qji  cos  (6  -  6  )  +  p,  cos  (6  -  6.)^  . 


We  have 


N 


u  cos  (6^ 


cos  (h*£ . 


6,) 


and  for  the  term  in  this  sum  we  first  average  over  the  electron, 
keeping  the  remaining  electrons  as  well  as  ions  fixed,  and  use  the  equation 
for  <cos(k»r  -  &)>  .  Next  we  average  over  the  remaining  N-1  electrons, 
still  keeping  the  ions  fixed,  and  add  the  identical  N  terms  j  =  1  to  N  . 
This  gives 


<p  cos  (6  -  6:)  > 

e  '  Pi 


[I<p 


cos  (6^ 
e  e 


6i)> 


P, 


=  -  d^d  +  a^)'Vi  .  (15) 

where  the  subscript  p.  indicates  that  iiiis  quantity  is  kept  fixed.  We 
finally  square  the  expression  in  Equation  (15),  average  over  the  positions 


jji - 

The  fallacy  in  Kahn's  derivation  appears  to  lie  in  the  use  made  of  his 
Equation  (17):  In  this  equation,  two  expressions  for  <pr>  occur,  one 
for  N  electrons  and  one  for  N  +  1  ,  which  should  be  taken  at  constant 
density  N/V  (not  at  constant  volume  V).  Since  second  differences 
occur  ihi  subsequent  equations,  the  use  of  a  constant  volume  V  is  not 
justified. 


-11- 


of  all  the  ions  as  well  and  use  the  explicit  expression  in. Equation  (13) 
for  <ii?>  to  obtain  the  desired  result 

=^HgCos(6^  -8i)>pVzNeV(l+a2)'tl  +  (Z+l)o2^-^  , 


e  =  <u^>  -  e  .  =  Ne^  (1  +d^)'' 

ee  '^e  ei 


If  the  ion  mass  M  is  large  compared  with  the  electron  mass  m  , 
the  expressions  in  Equation  (16)  have  the  following  physical  significance. 

The  ions  move  slowly,  so  that  varies  slowly  with  time,  and  the  fre¬ 
quency  Fourier  transform  of  p?  has  a  narrow  spread  as^  does  the  Fourier 

transform  of  the  square  of  that  part  of  p  which  remains  in  phase  with  p  . 

®  i 

The  full  Fourier  transform  of  Pg(t)  (which  we  shall  analyze  in  detail  in 
the  next  section)  thus  contains  one  part,  with  a  narrow  frequency  spread, 
whose  integrated  intensity  is  given  by  ©gj  •  The  remaining  part,  repre¬ 
senting  the  electron  density  fluctuations  which  are  not  correlated  with  the 
ions,  has  a  wide  frequency  spread  characteristic  of  electron  thermal  velocity 
Doppler  broadening  and  integrated  intensity  ^iven  by  ®g  g  •  For  a  <^<  1 

the  dominant  part  is  ©gg  and  the  integrated  intensity  of  the  narrow  part 

2  2  4 

IS  small;  Ogg  %  Ne  and  0^^  ~  ZNe  a  .  For  a»  1  the  dominant  part 

is  0  and 
e  1 


®ee  ' 

In  deriving  Equation  (10)  from  Equation  (8)  for  a  particular  wave 
mode  k  we  carried  only  the  leading  term  in  an  expansion  in  powers 
of  1//CT  and  thereby  omitted  all  terms  which  involve  any  other  wave 
mode  q  ^  k  ,  The  approximation  made  thereby  is  equivalent  to  the  so- 
called  "random  phase  approximation,  "  which  neglects  correlations  between 
different  wave  modes.  Some  deviations  from  this  approximation  could  also 
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be  calculated  with  methods  similar  to  those  of  the  present  section.  For 


instance,  with  p.,  ahd  p,  .  fixed,  the  average  of 

.  k,e- 


u  lx  cos  (6-  6  ^  -  6^,.) 

^^qt^q't  '  qt  q't' 

for  q  q'  =  k  coiild  bc[  evaluated.  Using  such  expressions  and  keeping 
terms  of  second  power  In  1//CT  in  Equations  (8)  and  (9),  corrections  to 


Equations  (10)  and  (13)  could  then  be  obtained.  For  a  «  1  the  leading' 
correction  to  Equation  (10)  is  probably  of  relative  order  e^k//CT  ~  (Aa)  ^  , 
the  correction  to  of  relative  order  o/A  ,  where  a  and  A  are 


defined  in  Equations  (4)  and  (5).  For  o  »  1  the  leading  correction  both 
to  Equation  (10)  and  to  is  probably. of  relative  order  (Ao^)  ^ 


III.  THE  FREQUENCY  SPECTRUM 

We  have  so  far  evaluated  only  the  root-mean  square  average  of  the 
quantity  Pg(t)  defined  in  Equation  (1)  and  now  wish  to  calculate  its  time 
dependence,  or  rather  Us  frequency  Fourier  transform.  In  this  section 
we  neglect  collisions  entirely  (mean  free  path  much  larger  than  both  k"^ 
and  D)  and  assume  that  the  only  forces  acting  are  those  of  the  electric 
field  E(r)  ,  given  in  Ecjuation  (7),  arising from’the  Charge  density  fluctuations 
themselves.  We  again  assume  that  A  »  1  and  also  that  nk  ^  »  1  .  In 
this  case  the  use  of  a  Boltzmann  equation  for  a  Boltzmann  distribution 
function  f(r  ,  V,  t)  is  in  general  justified.  Such  an  equation  was  used  by 

7 

Bhatnagar,  Gross,  and^Krook  and  the  spatial  and  frequency  Fourier  trans¬ 
forms  of  the  distribution  function  f  evaluated.  We  shall  use  a  method 
similar. to  theirs,  but  We  shall  have  to  take  account  of  the  discrete  nature 

of  the  electrons  more  explicitly  since  we  wish  to  retain  also  terms  in  an 
2 

expression  for  pg(t)  which  are  proportional  only  to  N  ,  the  number  of 
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2 

electrons  in  the  volume  V  ,  rather  than  N  . 

We  consider  a  fixed  value  of  the  wave  vector  k  ,  take  its  direction 
as  the  2  axis,  call  the  z  component  of  velocity  v  and  shall  omit  the  sub¬ 
script  k  in  ,  etc.  We  define  a  quantity  (t)  by 


d  V  O'  (t) 
ev 


(17) 


(v) 

where  2  denotes  summation  over  all  electrons  whose  velocity  lies 
between  v  and  v  +  dv  .  Our  quantity  o'gyit)  is  essentially  -e  times 
the  spatial  Fourier  transform  for  wave  vector  k  of  the  Boltzmann  dis¬ 
tribution  function  f(r,  v ,  t)  .  The  Boltzmann  equation  reads 


where  the  right-hand  side  represents  the  contribution  from  those  electrons 
whose  velocity  was  below  v  previously  but  passed  into  the  velocity 
regions  between  v  and  v  +  dv  .because  of  the  acceleration  by  the  electric 
field  E  (minus  those  that  have  passed  beyond  v  +  dv);  E  is  given  by 
Equation  (7).  We  assume  A=  nD^  »  1  ,  nk"^  »  1  ,  and  the  absence  of 
any  large-scale  macroscopic  deviations  from  the.rmal  equilibrium,  and 
replace  f  in  the  right-hand  side  of  the  Boltzmann  equation  by  the  equili¬ 
brium  distribution  function 


f(r  ,  V,  t)  —  nFp(v) 


1  2 

F^(v)  =  /2/CT 


(18) 


where  n  =  N/V  is  the  electron  particle  density.  This  gives 


B(t  /8t  +  ikvo- 
e  V  e  V 


i(4Tre^n//CT)(v/k)  p  F  (v)  , 

e 


(19) 


where  is  Ppj^^  for  our  fixed  value  of  k  .  This  substitution  is  equivalent 
to  the  "random  phase  approximation,  "  since  we  have  neglected  fluctuations 


in  f(£,  V  ,  t)  with  wave  vector  q  together  with  components  for  E  in 
Equation  (7)  with  wave  vector  k  -  q  . 

We  define  next  the  frequency  Fourier  transforms,  or  rather 
Laplace  transform^,  of  the  time -niependent  quantities  <rgy{t)  and  pg(t)  , 


dtv^^  .  Q^(co)=/  dtp^aie'^'^"^^^* 


where  Y  is  a  real,  positive,  infinitesimally  small  constant.  In  a  radar 
experiment,  where  tHe  frequency  spectrum  of  '|pg(t)|^  is  obtained  during 
a  large  but  finite  time  interval  (2Y)~^  ,  a  quantity  essentially  like 
|Qe(“)|^  is  measured  and  we  also  have  for  the  time  average  of  jpgp  > 


we  derive  from  Equation  (19)  the  relation 

=  (u  +  kv- W)'*C- »  (■4ireV/CT)(v/k)F^(v)Q^(uO  , 

(22) 

where  ,  p^,  q/^^  and  are  defined  by  equations  analogous  to 

Equations  (17)  and  (20)  but  with  all  the  charges,  electrons,  and  positive 
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ions  of  charge  Ze  ,  included.  For  the  ions  alone  one  obtains  an  equation 
sinnilar  to  Equation  (22)  with  ne  replaced  by  Zne^  and  with  a  distri¬ 
bution  function  F^(v)  occurring,  which  is  defined  as  in  Equation  (18) 
but  with  the  ion  mass  M  replacing  the  electron  mass  m  .  In  this  ex¬ 
pression  we  also  allow  the  ion  temperature  T.  to  differ  from  the  electron 
temperature  T  .  For  the  electron  and  ions  combined  we  then  find 


=  {«  +  kv-iY)-*[:-io-^y{0)-(4nneVkKT)(F^+ZTTrlFi)Qj(«)l  . 


Summing  Equation  (23)  over  all  velocity  groups  v  and  using  the 
fact  that  Q  ,  we  obtain  an  explicit  expression  for  Q^(cj)  in  terms 

of  the  quantities  the  initial  time  t  =  0  .  The  terms  involving 

Fg(v)  and  Fj(v)  are  smoothly  varying  functions  of  v  ,  and  we  can  replace 
the  summation  over  v  by  an  integration.  In  the  term  involving  (r^(0)  , 
however,  we  must  be  careful  to  preserve  the  discreteness  of  the  summation 
and  of  expressions  like  Equation  (17)  and  write  this  term  as  a  summation 

over  individual  electrons  and  ions.  This  gives 

{ N  N/Z 

Q^(w)  =  ieO  +  Gg(w)  +  G-(w)3]  (w  + kVj  -  iY)"^e  "^-Z  («  +  kvj  -  iy)  "^e 


where 


G  (w)  =  /  (4ffe^nv/k/CT)(a)  +  kv  -  i'Y)"^F  (v)dv  , 

6  /  C 


Gi^u)  =  /  (4tt  Ze^nv/k/CT.)(w  +  kv  -  •Y)"^F.(v)dv  , 


and  z.  ,  Z.  denote  the  position  of  the  electron  or  ion,  respectively,  at 
J  J 

time  t  =  0  (and  v.  is  the  corresponding  velocity).  After  summing 
J  • 
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Equation  (22)  over  velocity  v  and  making  use  of  Equation  (25),  we  also 
find 

N  -u  N/Z 

-1  k  ^  -i  K  z, 


Qg((i))  =  ie 


Y"  e  j 

1  +  Gg  +  Gj  /  w  +  kVj  -  iY  1  +Gg  +G-  /  w  +  kvj  -  iY 


j  =  1 


j  =  1 


(26) 


Equation  (26)  expresses  Qg(&)  explicitly  in  terms  of  the  positions 

and  velocities  of  all  the  electrons  and  ions  at  the  initial  time  t  =  0  .  If  we 

were  to  put  Z  =  0  and  replace  the  summation  in  Equation  (26)  by  an  integral 

over  (T  (0)dv  ,  we  essentially  wo>uld  obtain  Equation  (44)  of  Bhatnagar 
ev 

7 

et  al.  ,  and  we  shall  return  to _sudh  am  equation  in  Section  V.  At  the  moment, 
however,  we  want  an  expression  for  jQ^(w)p  under  conditions  of  thermal 
equilibrium  (except  that  Tj  may  differ  from  T)  and  therefore  take  the 
modulus  squared  of  the  right-hand  side  of  Equation  (26)  and  average  over 
initial  conditions.  In  the  double  summation  over  particles  j  and  I  we 
separate  out  the  terms  with  j  =  /  and  obtain 
N 

V  - rrrl  )=  N 


-i  k  z.. 
J 


j  =  i 


cj  +  kv.  -  iY 
J 


00 


(27) 


The  quantity  expressed  as  <  >  V2  i^^dicates  an  average  over  the 

positions  of  the  two  electrons  with  their  velocities  kept  constant  and  is  of 
the  same  order  of  magnitude  as  the  expression  in  Equation  (14),  With  the 
density  fixed,  the  second  term  on  the  right-hand  side  of  Equation  (27)  is 
then  proportional  to  N  ,  just  as  the  firs^  term,  but  as  Y  0  the  second 
term  tends  to  a  constant  limit  whereas  the  first  term  is  proportional  to 
Y"^  .  Keeping  only  terms  of  order  Y'^  (see  Appendix  I),  we  can  neglect 
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the  second  term  in  Equation  (27)  and  evaluate  the  integral  in  the  first 
term  and  obtain  from  Equation  (26) 


(v//7rNe^<|Qg(a))|2>=  |l  +  +  G.|'^k-lyi  +  qp  (-u,A)  +  Z  |g^|^  F.  (-u)A)  J  . 

(28) 

where  G  is  defined  in  Equation  (25)  and  F  in  Equation  (18). 


IV.  RESULTS 


Equation  (28)  is  the  essential  result  in  its  most  general  form.  For 

Q 

Z  =  T/T.  =  1  this  result  has  also  been  obtained  by  Dougherty  and  Farley 
and,  for  a  »  1  ,  by  Fejer.  Using  Equations  (54)  and  (56)  of  P.  L.  Bhatnagar, 
E.^P.  Gross  a:nd  M.  Krook,  the  expressions  for  G  in  Equation  (25)  can  be 
expressed  in  terms  of  tabulated  functions  (see  Appendix  II)  : 


CM  -= 


1  -  £{xy~‘islrr  xe 


-X 


=  w/a)^  .  =  (2k^A:T/m)i  ,  (29) 


G.(w)  -  =  (ZTa  /T.)  |l  -f(y)  -  i^fnye'y 
2 

f(x)  =  2 xe“^  /  e*”  dt 


21 


,  y  =  u)/u.)  ,  Uj  =  , 

(30) 


Using  tables  of  the  integral  in  Equation  (30)  we  have  tabulated  f(x)  .  These 
values  are  given  in  Table  I.  For  x  <  1  the  Taylor  series  (convergent  for 
all  X  )  converges  rapidly  , 


f(x)  =  2x^  {1  -  (2/3)x2  +  (4/l5)x^+  ...  -  (2x^)*'  [3-5--(2n+  1)] 
For  X  »  1  we  have  the  asymptotic  expansion 


+  ...} 
(31) 
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2  ^ 

Table  I.  f(x)  =  2xe"^  J 

o 

2 

dt . 

X 

f(x) 

X 

f(x) 

x 

f(x) 

.  05 

.  0050 

1.  30 

1.2568 

3.  10 

1. 06392 

.  10 

.0199 

1.  35 

1.2696 

3.20 

1. 05895 

.  15 

.  0443 

1.40 

1.2782 

3.  30 

1. 05461 

.20 

.  0779 

1.45 

1.2831 

3.  40 

1. 05075 

.25 

.  1199 

1. 50 

1.2848 

3.  50 

1. 04735 

.  30 

.  1696 

1. 55 

1.2835 

3.  60 

1. 04429 

.  35 

.  2259 

1.60 

1.2797 

3.  70 

1. 04157 

.40 

.  2880 

1.65 

1.2741 

3.  80 

1. 03909 

.45 

.  3545 

1. 70 

1.2668 

3.  90 

1. 03682 

.  50 

.  4244 

1.  75 

1.  2580 

4.  00 

1. 03478 

.55 

.  4966 

1. 80 

1.  2483 

4.  10 

1.03293 

.60 

.  5697 

1. 85 

1.2380 

4.20 

1.03119 

.  65 

.  6428 

1.  90 

1.  2272 

4.  30 

1. 02965 

.  70 

.  7147 

1. 95 

1. 2160 

4.40 

1.02815 

.  75 

.  7845 

2.  00 

1.2057 

4.  50 

1.02679 

.  80 

.  8514 

2.  10 

1.  1844 

4.  60 

1. 02554 

.  85 

.  9145 

2.20 

1.  1639 

4.  70 

1. 02437 

.  90 

.  9733 

2.  30 

1.  1457 

4.  80 

1.  02328 

.95 

1.  0273 

2.  40 

1.  1295 

4.  90 

1. 02229 

1. 00 

1.  0762 

2.  50 

1.  1154 

5.  00 

1. 02134 

1.05 

1.  1196 

2.  60 

1.  1032 

1.  10 

1. 1577 

2.  '70 

1. 09281 

1.  15 

1.  1902 

2.  80 

1. 08389 

1. 20 

1.2175 

2.  90 

1.07622 

1.25 

1.2396 

3.  00 

1. 06959 
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f(x)  .  1^  (2x^)‘^  Cl  +  3/(2x2)  +  15/(4  x^)  +  .  .  .  C 


(32) 


For  intermediate  values  of  x  ,  the  function  f  (x)  is  plotted  in  Figure  1  , 
as  is  the  function  obtained  from  the  first  two  terms  on  the  right-hand  side 
of  Equation  (32)  (dashed  curve).  The  velocity  distribution  functions  F  in 
Equation  (28)  can  be  written  explicitly  as 


k"^F  =  e^/ ^/Tr 

e  '  e 


and 


k"  F.  =  e  ^  / 'n/ttu. 


The  constants  and  represent  Doppler  spread  frequencies 
characteristic  of  thermal  velocities  of  the  electrons  and  ions,  respectively. 
For  most  cases  of  practical  interest  m  «  M  and  T- 51  T  so  that 

Ti  =  w./cj  =  (mT./MT)i  «  1  .  (33) 


In  this  case  a  good  approximation  (except  for  some  special  cases  discussed 
below)  can  be  given  for  Equation  (28),  in  terms  of  a  single -parameter 
family  of  functions  F  of  one  variable,  as  follows.  The  first  term  in 
Equation  (28)  involves  F-  OC  e"^  and  is  of  most  interest  for  |x|  oe  1  . 


Disregarding  the  narrow  region  |x|  =  |yp  ,  we  have  |y|  »1  for 

this  term  and  ^-ZTa^/2T^y^  can  be  neglected  compared  with  unity 

and  with  .  The  second  term  in  Equation  (28)  involves  F^  oce'^  and 

is  unimportant  if  jy  =  jxj  •  In  the  important  regions  we  thenhave 

jx  <<  1  and  Gg  J^a^(l  +  i^/i^x)  .  Neglecting  also  the  term  isTnx  ,  we  obtain 
.* 

2 


finally 


(34) 


T.(l  +  0^-) 


^This  result  has  been  stated  previously  by  E.  E.  Salpeter. 
Equation  (1)  l/2  it  should  read  l/2‘trp^0^  . 


10 


In  his 


.-20- 


1.30 


rjx)  =  e 


f(x)j  ^  +  tr 


4  2  -Zx 
axe 


where  a  is  defined  in  Equation  (5)  and  x,  y  in  Equation  (29) •  Values 

for  Fjj  are  given  in  Table  II  for  several  values  of  a  . 

Each  function  F^fx)  is  even  in  x  .  It  is  plotted  for  positive  x  in 

Figure  2  for  a  =  0  ,  0.5  ,  1,  2,  3  and  4  .  For  o  «  1  the  function  is 

2 

close  to  the  Gaussian  F  (x)  =  e  .  This  is,  of  course,  the  characteristic 

D 

Doppler  spread  spectrum  for  noninteracting  electrons.  For  a'»  1  ,  on  the 
other  hand,  F^  (x)  has  a  very  sharp  maximum  near  x  =  ±.  x^  ,  where  Xq 
is  the  solution  of  the  dispersion  relation 


f(x  )-l  =  a 
o 


For  X  very  near  x^  ,  F^^  (x)  can  then  be  approximated  by  the  Lorentzian 


shape 


r^(x)  %  ia^e"*o 


r  ?  A  221-1 

4{x  -  Xq)  +  (InTu  a^e'^°)  J 


for  a  »  1  ,  where  we  have  used  the  approximate  relation  f  -  1 
in  evaluating  coefficients.  If  Equation  (36)  is  solved*  approximately  by 
using  the  first  two  terms  in  the  asymptotic  expansion,  Equation  (32),  we 


obtain 


x2%|(a2  +  3)  .  0)^  =  +  3/CTkV' 


This  expression  for  is  the  well-knownll  dispersion  relation  for  longi- 

o 

tudinal  (electrostatic)  plasma  oscillations.  The  Lorentzian  shape  of 

♦Equation  (36)  also  has  a  second  solution  with  Xq  ~  1  .  This  solution  is  of 

no  interest  since  F’  ~  «  1  . 

a  - 

♦♦The  relevance  of  such  plasma  oscillations  to  the  radar  backscatter  problem 
was  first  pointed  out  by  Akhiezel?  et  al,  , T .  who  used  ;a  model  in  which  the 
Lons  are  replaced  by  a  uniform  charge  distribution. 
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Figure  2.  The  function  r^(x)  plotted  against  x  for  u  =  0,  0.5,  1,  2.  3  and  4 .  The  vertical  lines  near  the  top  of  the  Figu 
denote  x^  ,  Equation  (38),  for  a  =  2,  3  and  4  .  The  dashed  curve  denotes  r**^x).  Equation  (4f),  for  0=1. 


Equation  (37)  is  characteristic  of  the  resonance  spectrum  for  a  long- 
lived  oscillation.  The  width  of  the  spectrum,  the  expression  in  the  second 
parenthesis  in...  Equation  (37)  comes  from  the  so-called  Landau  (or 
"drift")  damping,  which  is  contributed  by  those  few  electrons  in  the  tail 
of  the  Maxwell  distribution  whose  velocity  equals  the  (very  large)  phase 
velocity  of  the  plasma  oscillation.  As  a  increases,  the  width  decreases, 
and  the  maximum  of  r^(x)  increases  sharply,  even  though  the  integrated 
intensity, 

00 

n  J"  r^(x)  dx  %  a  ^  ,  (39) 

-00 

decreases.  It  should  be  remembered  that,  for  a  practical  problem, 
collisions  also  contribute  a  very  small  width  to  the  spectrum,  which 
dominates  the  Landau  damping  for  very  large  values  of  a  and  that  small 
slow  variations  of  the  over -all  electron  density  will  vary  cUp  and  broaden 
the  spectrum.  As  the  discussion  in  Section  II  shows,  the  integrated  in¬ 
tensity  in  Equation  (39)  should  not  depend  on  collisions  or  on  the  width. 

We  have  discussed  so  far  only  the  first  term,  involving  r^(x)  , 
in  Equation  (34)  which  represents  the  part  of  the  frequency  spectrum  that 
is  important  at  large  frequencies,  of  the  order  of  or  of  co^  ,  and 
whose  integrated  intensity  is  given  by  in  Equation  (16).  We  turn 

now  to  the  second  term,  involving  rp(y)  ,  in  Equation  (34)  which  is  im¬ 
portant  only  for  small  frequencies  and  whose  integrated  intensity  is  given 

Equation  (l6).  For  a  «  1  we  have  %  Zo^  ,  0  «  0 

61  61  6666 

but  the  width  of  the  second  term  is  smaller  by  a  factor  of  n  =  w./wg  «  1 
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than  that  of  the  first,  and  the  peak  of  the  second  term  will  dominate  for  small 
w  as  long  as  Z  »  ti  .  For  a  »  1  the  integrated  intensity  0^^  of  the 
second  term  dominates  that  of  the  first  term.  For  the  case  of  greatest 
interest,  Z  =  Tj  /  T  =  1  ,  we  then  have  P  =  1  as  a  -►  ob  and  (y) 
has  the  almost  flat-topped  shape  plotted  in  Figure  2.  In  this  case  we  have 
almost  complete  charge  neutrality;  the  electron  density  mainly  follows 
that  of  the  ions  which  can  change  only  slowly  and  leads  to  a  narrow  fre¬ 
quency  width  of  order  .  The  shape  of  Fj^y)  differs  from  the  Gaussian 
for  noninteracting  ions  because  electrostatic  potentials  of  order  /C  T 
are  set  up  by  the  requirement  that  the  electrons  follow  the  charge  density 
of  the  (slow)  ions. 

If  the  ion  temperature  is  lower  than  the  electron  temperature 
T  ,  as  well  as  a  »  1  ,  we  have  p  =  V  Z  T/T .  >  1  ,  and  the  "ion  com¬ 
ponent"  r'p(y)  also  has  a  Lorentzian  shape  like  Equation  (37).  This 
sharp  "resonance  curve"  represents  the  so-called  positive-ion  oscillations^ 
whose  frequency  is  the  same  as  that  of  a  plasma  osciUation  for  fictitious 
particles  with  the  ion  charge  and  mass  but  with  the  electron  temperature. 

If  the  quantity  q  ,  defined  in  Equation  (33)  ,  is  negligibly  small  and  there 
are  no  collisions,  the  width  of  the  frequency  spectrum  for  the  positive-ion 
oscillation  is  given  by  Fp(y)  ,  no  matter  how  large  p  is  .  However,  if 
(tiT^/ZT)  »  e“^  /2  ,  then  the  replacement  in  the  derivation  of  Equation  (34) 
of  (1  +  isTirx)  by  -a  is  not  justified,  and  the  actual  width, 

although  still  small,  is  larger  than  that  given  by  Fp(y)  . 

To  summarize  the  results  so  far,  for  the  most  important  case  of 

4 

T^  =  T  ,  Z  =  1  and  m  «  M  .  We  have  defined  the  dimensionless  para¬ 
meter  a  in  Equation  (5)  and  have  p  =  a(l  +  a^)~  ^  ,  t|  =  w./w  =  V m/M  . 

1  ©  » 
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The  frequency  distribution  is  given  by  Equation  (34)  with  x  and  y  defined 
in  Equation  (29).  The  integral  over  du  of  the  first  term  in  Equation  (34) 
is  ^/T^(l  +  a^)  ;  that  of  the  second  term  is 

a^(l  +  o^)  ^  (1  +  Zaf  )  ^  . 

The  sum  of  the  two  integrals  is 

<^^(1  +  a^)  (1  +  2a^)‘^ 

which  decreases  only  by  a  factor  of  two  as  a  goes  from  zero  to  infinity. 
The  function  r^(x)  is  even  in  x  and  is  plotted  for  positive  x  in 
Figure  2,  has  Gaussian  shape  for  a  =  0  ,  is  almost  flat-topped  for  a  =  1  , 
and  has  a  maximum  at  a  nonzero  value  of  x  for  larger  values  of  a  . 

j 

For  a  >  4  the  function  ^^^(x)  has  the  Lorentzian  shape  of  Equation  (37) 
with  a  very  sharp  peak  of  height  greater  than  unity,  and  Equation  (38)  is 
a  very  good  approximation.  For  a  <  (m/M)  '  ,  the  maximum  of  the 

first  term  in  Equation  (34)  is  larger  than  that  of  the  second  term.  This 
is  again  the  case  for  large  values  of  a'  (a  ~  5  for  M/m  ~  lO"^)  ,  but  it 

I 

should  be  remembered  that  in  a  practical  problem  such  as  the  ionospheric 
application,  there  are  other  causes,  besides  Landau  damping,  broadening 
the  "resonance  peak"  ,  and  the  actual  maximum  will  be  lower  than  that 
given  by  r^(x)  . 

y..  SOME  DEVIATIONS  FROM  THERMAL  EQUILIBRIUM 

We  have  discussed  so  far  only  cases  in  which  complete  thermo¬ 
dynamic  equilibrium  holds  (except  that  the  ion  and  electron  temperatures 
Tj  and  T  may  differ).  We  now  consider  one  very  special  kind  of  small 
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deviation  from  equilibrium.  We  assume  that  equilibrium  has  been  established, 
but  that  at  some  time  t  =  0  ,  the  electron  and  ion  charge  densities  are  sud¬ 
denly  both  altered  in  a  nonhomogenous  manner  (although  the  charge  is  kept 
neutral)  by  some  external  agent.  This  might  be  accomplished,  for  instance, 
by  the  sudden  passage  of  fast  ionizing  particles  with  a  patchy  spatial  distri¬ 
bution.  The  newly  created  patchy  electron  charge  distribution  is  assumed  to 
be  small  compared  with  the  uniform  density  n  ,  but  its  spatial  Fourier  trans¬ 
form  (0)  =  -  (0)  of  Equation  (1),  is  assumed  to  be  larger  than  p^^^ 

for  purely  thermal  density  fluctuations.  We  further  assume  the  absence  of 
collisions  and  A  »  1  in  Equation  (4)  ,  will  use  the  random  phase  approxi¬ 
mation  and  consider  only  one  particular  value  of  the  wave  vector  k  (and  drop 
the  subscript  k).*  We  assume  next  that  (0)  and  (t.  (0)  of  Equation 

(17),  are  smoothly  varying  functions  of  the  velocity  v  .  We  also  assume  that 
after  the  initial  time  t  =  0,  there  are  no  external  forces  or  disturbances 
(except  for  the  possibility  of  another  sudden  burst  of  ionization  after  the  effects 
of  the  original  disturbance  have  died  down). 

If  we  assume  that  the  disturbance  at  t  =  0  occurs  instantaneously, 
then  the  frequency  Fourier  transform  of  Pg(t)  is  again  given  by  Equation 
(26),  but  we  can  replace  the  summations  by  integrations  over  the  smoothly 
varying  functions  (0)  and  (0)  .  We  specialize  further  by  assximing 

that  (0)  is  proportional  to  the  Maxwell  distribution  function  for  ions 
at  temperature  T  =  T-  and  that  (r  (0)  is  proportional  to  the  Maxwell 

1  V  c 

distribution  for  electrons  at  temperature  T6"^  where  6  i  1  .  Carrying 

^Because  of  the  absence  of  collisions  our  case  is  quite  different  from  those 
involving  turbulence  (with  a  short  mean  free  path)  where  large  eddies  feed 
small  ones  and  also  depends  on  ,  pjQg(O)  with  q?^k;  see,  for 

instance,  Villars  and  Weisskopf^®  and  Silverman^^. 
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out  the  integrations  over  dv  we  find 


QgN  = 


^  +  °i  f.  .  £(x6)  1  .  °e  r.  -  /  ,  f(y)  1 

—  -  le  -f  -■  r-  + -  le  +  -7: — 

i  “g  x6\MtJ  L 


where  a  is  a  constant.  Using  the  fact  that  m  «  M  we  can  simpli|^y  this 
expression  as  we  did  in  Section  IV  to  obtain 


Q^Mp  ocr^^^(x)a);^  +  Zd^l+ 

r^i^(x)  =  e"^^  +  e*^f^(x)/iTX^  r^(x)  .  (41) 

As  Equation  (41)  shows,  the  function  r^^^(x)  decreases  much  less 

®  2 
2  —X 

rapidly  for  large  x  than  the  function  r^(x)  ,  as  x"  rather  than  as  e 
This  slow  fall -off  is  due  to  our  special  assumption  of  a  sudden  onset  of  the 
disturbance  which  contributes  Fourier  components  of  large  frequency.  If 
the  onset  occupies  a  finite  time  duration  T  ,  as  it  woujd  in  practice,  our 
Equation  (41)  breaks  down  for  «>  T"^  and  the  actual  spectral  intensity 
would  be  lower  than  in  our  approximation.  Note  also  that  |Og(«)p  in 
Equation  (41)  is  independent  of  ,  the  length  of  time  over  which  the  fre¬ 
quency  spectrum  is  accumulated  (as  y  0),  rather  than  being  proportional 
to  y”^  ,  as  is  the  expression  in  Equation  (34)  for  the  case  of  thermodynamic 
equilibrium.  Since  we  assumed  only  the  creation  of  a  single  external  dis¬ 
turbance,  which  dies  down  in  a  finite  time  period,  the  time  periods  beyond 
that  point  do  not  contribute  to  Equation  (41). 

For  Q  1  ,  (and  6  1  ),  the  peak  values  of  do  not  differ  very 

greatly  from  those  of  .  For  a  »  1  the  second  term  in  Equation  (41) 

contains  F^P(y)  .  The  dashed  curve  in  Figure  2  depicts  this  function  F^p(x) 
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which  is  seen  to  be  similar  to  rj(jc)  except  for  its  longer  tail.  The  first 
term  in  Equation  (41)  ,  however,  behaves  rather  differently  for 

a  »  1  (even  with  6  ~  1)  .  It  has  the  Lorentzian  shape  of  the  expression 

in  square  brackets  in  Equation  (37)  but  the  very  small  multiplying  factor 

2  { 1 ) 

x^e  °  is  missing.  The  integral  of  over  x  thus  increases  with  a 

-ZA  _2 

roughly  as  e®  rather  than  decreasing  as  a  .  Physically  this  means 

that  our  assumed  sudden  onset  of  the  disturbance  can  excite  a  plasma 

oscillation  no  matter  how  much  larger  cjp  is  than  ,  but  this  oscillation 

persists  for  a  length  of  time  (the  inverse  of  the  Landau  damping  frequency 

2/2 

width)  which  increases  with  a  as  e*^  '  .  It  should  be  noted  again  that, 

in  practice,  collisions  will  put  an  upper  limit  to  the  persistence  time  of 
the  plasma  oscillations  and  that  the  excitation  of  the  oscillations  would  be 
strongly  depressed  if  the  onset  time  T  of  the  disturbance  is  large  compared 
with  the  oscillation  period  Wp^  .  The  assumptions  in  this  section  were 
chosen  not  so  much  because  they  are  physically  reasonable  but  because 
their  consequences  follow  readily  from  previous  work  in  this  paper. 
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APPENDIX  I. 


Evaluation  of 


00 


F(v)  dv 

(w  +  kv)^  +  Yp 


Setting  X  =  w  +  kv  we  get 


I(Y).  = 


00 


F(^) 

+  Y^ 


oo 


dx 

k 


Now  let  X  =  Y  tan  0  ,  then 


Aften  integration  we  get 


Note  that  for  any  sufficiently  well-behaved,  function  F(v)  ,  —  (-w/k) 
is  the  only  term  in  I(v)  which  diverges  as  Y  -  0  . 
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APPENDIX  II.  Proof  of  Equation  (29) 


00 


G^(w)  = 
e 


2, 

a  kv 


Oj^kv-iY  ® 


F^(v)!dv 


-00 


Setting 


o)  =  v/Vz/Ct/  m 


w  m  _  (0 

k  «2/CT  ■  uth 


we  get 


00 


°e<“>  =  W 


00 


u  e  du 


X  +  u  -  iy 


Vir 


ue~^  du 
u  -  0 


-00 


-00 


(1) 


where  0  =  iY  -  x  .  Since  the  integrand  is  analytic  except  at  u  =  0  we  can 
change  the  path  of  integration  to  the  contour  c  shown  below. 


Hence  Gg(co)  is  resolved  into  a  contribution  from  the  small  semi¬ 
plus  the  principal  value  of  the  integral  along  c  . 

To  calculate  the  contribution  (Ij)  to  from  the  semi 

integration,  we  set  u  =  0  +  r  e^*^  .  Then 
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The  principal  value  contribution  (I2)  yields 
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Proof  by  identity 


Now  as  Y  -*  0  ,  0  -♦-x: 


2  r  2  .^2 

G^M  =  Ij  +  D  "  2xe'^  /  e  6t  -  i  xe  1 


=  a^Ql  >  f(x)  -  i^/7^xe■^  '2 


This  is  just  Equation  (29). 
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